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Abstract. The binding of polarons, or its absence, is an old and subtle topic. 
After defining the model we state some recent theorems of ours. First, the transition 
from many-body collapse to the existence of a thermodynamic limit for N polarons 
occurs precisely at U = 2a, where U is the electronic Coulomb repulsion and a is the 
polaron coupling constant. Second, if U is large enough, there is no multi-polaron 
binding of any kind. We also discuss the Pekar-Tomasevich approximation to the 
ground state energy, which is valid for large a. Finally, we derive exact results, not 
reported before, about the one-dimensional toy model introduced by E. P. Gross. 



1. Definition and previous rigorous results 

The large polaron, first considered by H. Frohlich [6] in 1937, is a model of an 
electron moving in three dimensions and interacting with the quantized optical modes 
of a polar crystal. (It is called 'large' because the size of the electronic wave function 
is large compared to the crystal lattice spacing, so a continuum approximation is 
appropriate.) It is also a simple quantum field theory model and over the years has 
been used as a testing ground for various approximations. 

In suitable units, its Hamiltonian is 

#W=p 2 +/ a \k)a(k)dk+^- [ -*-[a(k)exp(ik-x) + h.c.]dk. (1.1) 

Jm 3 v 2 7r Jr3 |k| 

This Hamiltonian acts in the Hilbert space L 2 (IR 3 ) <g> J 7 , where J 7 is the bosonic Fock 
space for the longitudinal optical modes of the crystal, with scalar creation and anni- 
hilation operators tv(k) and a(k) satisfying [a(k), a'(k')] = 5(k — k'). The momentum 
of an electron is p = — z'V, and the coupling constant is a > 0. (Other authors have 
used a different convention, where a is replaced by a/y/2 [6j [I].) 

The ground state energy E^(a) is the infimum of the spectrum of H^ 1 '. Because 
of translation invariance, E^'(a) cannot be expected to be an eigenvalue, and indeed 
it is not jU [7]. The following rigorous results concerning E^(a) will be important in 
our analysis. 

(i) For all a, 



-a — a 2 



/3 < E {1 \a) < -a. 
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These upper and lower bounds are in [HI El [To] and [23] , respectively. As a 
consequence, E ( - 1 \a) ~ —a for a small, 
(ii) Pekar [25], using a product function, showed that 

E {1) (a) < -C P a 2 , 

for all a. The lower bound 

£ (1) (") > -C P a 2 - const a 9/5 

for large a was proved in [21]. It was proved earlier in [2], but without the a 9//5 - 



error estimate. Here, Cp = 0.109 [24] is the number determined by Pekar's 
variational problem for the electron density, 

C p = mU |V^x-/Y ^W |2| ^ |2 dxdy: M2 = l}. (1.2) 
Uir3 J,/k3 xIR 3 |x-y| J 

The minimizing ip is unique up to translations (and a trivial phase) |16j . 
(hi) There is a representation for E^^a) in terms of path integrals. In terms 



a) 



of the partition function Z^ l \T) = Trexp( — TH^}, one has £ (1 
— limr-i-oo T^ 1 log Z^(T). (Strictly speaking, Z^(T) does not exist because 
of the translation invariance of Hq and the infinite number of phonon modes. 
These technicalities can be handled by inserting appropriate cutoffs, to be 
removed at the end of the calculation [261 EI]-) It was shown in [3] that af- 
ter one integrates out the phonon variables, Z^'(T) has a functional integral 
representation 

Z«\T) = /^exp [g n T fl^±- ] , (1.3) 
y [2 7 Jo |x(i) -x(s)|_ 

where <i/z is Wiener measure on all T-periodic paths x(t). (In physics notation 
dfi = exp(— J Q T x(t) 2 dt) (ipath. Strictly speaking, t — s has to be understood 
modulo T, but this is irrelevant as T — > oo.) 

2. The multi-polaron problem and new results 
Of great physical interest is the binding energy of N polarons, with Hamiltonian 

N „ N 

H { v ] = VV + / a f (k)a(k) dk + Y] / 7^[a(k) exp(zk • Xj ) + /i.e.] dk 
i= i ■/ V2vr j=1> / |k| 

+ u J2 i^-^-r 1 (2.i) 

i<i<j<N 

and ground state energy E^\a). We ignore Fermi statistics for the electrons, because 
its imposition changes things only quantitatively, not qualitatively. The Coulomb 
repulsion parameter U > is equal to e 2 . The derivation of Hy in [6] implies that 
U > 2a, and this is crucial for thermodynamic stability, as we shall see. 
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The generalization of ( 11. 3ft is 



r T r T e-^dtds r T dt 



Jo Jo \Mt) - X£ J° l X ^^ ~ Xj '^l 

"(2.2) 

where d//^ is Wiener measure on all T-periodic paths (xi(t), . . . , xjv(t)). This is 
relevant for us since E^(a) = — lim T _ s>00 T _1 log Z^\t). 

The generalization of the Pekar approximation (II. 2p to the iV-electron case is the 
minimization of the following Pekar- Tomasevich functional for normalized functions 

tfj on M. 3N , 

f I W^dX + U^f mX)l2 dX-a[[ e*VLPjW dxdy, 

~[ JR3N J R 3N |Xj - Xj| J JR3xR3 l X _ y| 



where dX = Yik=i ^ Xfc - The density p^ of ip is defined as usual by 

- f — 
p v ,(x) = '^ / • • • ,x, . . . , xjv)| 2 • • -dxi ■ ■ -dx N 

with x at the i-th position, and c?Xj meaning that c?Xj has to be omitted in the prod- 
uct riA:=i^ x fc- This minimization problem is obtained from the original problem of 
minimizing (ty, H^^f) by restricting the allowed t/Vs to be of the form ip £g> $, where 
ip G L 2 (lR 3Ar ), $ G J 7 , and both ^ and $ are normalized. Since the Pekar- Tomasevich 
functional is the result of a variational calculation, its energy gives an upper bound 
to the ground state energy E^(a). 

2.1. Binding of multi-polaron systems. We first consider the bipolaron binding 
energy A £77 (a) = 2E^(a) — E^\a). For some time this was thought to be zero for all 
U > 2a, on the basis of an inadequate variational calculation, but it is now known [TJ 
to be positive for some U > 2a. The first question we address is whether AEjj(a) = 
for U sufficiently large. It is understood that the effective interaction induced by 
the phonon field for two polarons at large distances d is approximately Coulomb-like 
—2a/d, but this alone does not preclude binding. (The reason for 2a ■ distance -1 can 
be seen from the Wiener integral representation ( 12. 2ft , where there is a factor a/2, but 
the pair appears twice, and the integral f R e~^~ s ^ ds = 2.) The known existence 
of bipolarons for some U > 2a is an effect of correlations. It is a priori conceivable 
that correlations lead to an effective attraction that is stronger than Coulomb at large 
distances. If it were, for example, equal to {2a/ d) log(log(log(c?))), then this minuscule 
perturbation of Coulomb's law, which would be virtually undetectable by a variational 
calculation, would result in binding for all U . The absence of binding is a problem 
that has resisted a definitive resolution for many years. The following was proved in 
BE]: 
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Theorem 2.1 (Absence of binding for bipolarons). Let N = 2. For some con- 
stant C < 26.6, 

E ( J\a) = 2E {l \a) (2.3) 

whenever U > 2Ca. 

The constant 26.6 vastly exceeds the current, non-rigorous estimates of about 1.15 
[33| [27], so it is an open problem to find a more accurate rigorous bound. 

The existence of a critical repulsion strength for a bipolaron is consistent with the 
idea that the attractive interaction induced by the field is Coulomb-like, and therefore 
one expects that there is an iV-independent U c (a) such that there is no binding of any 
kind when U > U c (a). This was proved in jU |5] as well. 

Theorem 2.2 (Absence of binding for N polarons). For given a > there is a 
finite U c (a) > 2a such that 

E { u\a) = NE {l) (a) for all N > 2 (2.4) 

whenever U > U c (a). 



Remark 2.3. If U > U c (a), then given (I2.4p and any normalized i/j 

|^| V) > NE^(a) + (U- U c (a)) ^| J2 l x * ~ ^l" 1 ^) • (2.5) 

This inequality gives a quantitative estimate of the energy penalty needed to bring 
two or more particles within a finite distance of each other. In particular, it implies 
that for U > U c (a) there cannot be a normalizable ground state, even in a fixed 
momentum sector. Inequality (12.51) is not only true for our bound on U c (a), but also 
for the (unknown) exact value of the critical repulsion parameter. 

We state Theorem 12.11 separately for two reasons: One is that the proof is easier 
than for the general N case. The second is that we have an upper bound on U c (a) 
that is linear in a. While our iV-polaron bound is linear in a for large a, we have not 
achieved this linear bound for small a and this remains an open problem. 

2.2. Thermodynamic stability of multi-polaron systems. The second problem 
we consider is the existence of the thermodynamic limit. For large N, physical intuition 
suggests that E^\a) ~ —const N. This supposition is known to be false if U < 2a. 
Indeed, it was shown in [9] that, even with the Pauli principle, E^\a) ~ —const 
when U < 2a. Absent the Pauli principle, E^f\a) would behave even worse, as 
-const N 3 . It is also known [9] that E^\a) > -const iV 2 if U > 2a. The latter 
bound ought to be —const iV instead, and in [U [5] we proved that this is indeed the 
case for all U > 2a. 

Theorem 2.4 (Thermodynamic stability for U > 2a). For given U > 2a > 0, 
N~ 1 E i ff \a) is bounded independently of N. 
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Our lower bound on iV 1 E\ J (a) goes to — oo as U \ 2a, but we are not claiming 
that this reflects the true state of affairs. Whether limAr_>oo N~ 1 E^ l \a) is finite or 
not remains an open problem. There are partial results in the Pekar-Tomasevich 
approximation [9]. 

The linear lower bound from Theorem I2.4[ together with the sub-additivity of the 
energy [9], [2Ql Sec. 14.2], i.e., 

4" +M V)<4"V) + 4 M) («), (2.6) 

implies: 

Corollary 2.5 (Thermodynamic limit for U > 2a). For given U > 2a > 0, 
liniTv-j.oo N~ 1 E l {f\a) exists. 

For U in the range 2a < U < U c (a), there are bound states of an undetermined 
nature. Does the system become a gas of bipolarons, or does it coalesce into a true N- 
particle bound state? If the latter, does this state exhibit a periodic structure, thereby 
forming a super-crystal on top of the underlying lattice of atoms? This is perhaps the 
physically most interesting open problem. While particle statistics does not play any 
role for our main results, the answer to this question will crucially depend on particle 
statistics (Bose or Fermi) [30, 29J. 

3. Absence of bipolaron binding 

In order to give the flavor of our methods, we sketch the proof of Theorem 12.11 
as given in j4], ©Amer. Phys. Soc. The proofs of the other two theorems are also 
sketched in [4]. 

The proof of Theorem 1 is conveniently structured in 4 steps. 

Step 1. Partition of the interparticle distance: We fix a length £, whose value will 
later be chosen proportional to a" 1 , and partition the relative distance r = |xi — x 2 | 
between the particles into spherical shell-like regions of radial size 2 k ~ 1 £ < r < 2 k £ 
with k — 1, 2, . . .. This partitioning is one of the key points of our analysis. In addition 
there is the k = region, where the particle separation is between zero and £. Because 
of the uncertainty principle these regions have to overlap a bit, but this can be easily 
handled, and we ignore it for the sake of simplicity. There is a kinetic energy cost for 
localizing the particles according to this partition, which is c\2~ 2k £~ 2 in the shell k. In 
the next step we look at the energy of the particles localized to one of these shell-like 
regions. 

Step 2. Further localization for well- separated particles: For k > 1 we further 
localize the particles into individual boxes of size 2 k ~ 3 £. This costs another localization 
error C22~ 2k £~ 2 . Because the separation exceeds 2 k ~ l £, the two particles cannot be in 
the same or neighboring boxes. From the path integral (12. 2p . but now with the Xj(t)'s 
constrained to their respective boxes, we see that the separated particles feel an effective 
Coulomb-like attractive potential. However, this can contribute at worst — c 3 a2~ k £~ 1 
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to the energy. But the Coulomb repulsion is at least U2 k l 1 , which implies that the 
total energy exceeds 2E {1) if 

U2~ k r 1 > c z o2~ k f x + (ci + c 2 )2" 2fc r 2 . 

If this inequality holds for k — 1, it holds for all k > 2 as well. Thus, if we can deal 
with the k = region, we will establish that binding is not possible if 

Ua' 1 > c 3 + (ci + c 2 )/(2£a) . (3.1) 

5*tep 5. TTie region of no minimal separation: In the k = region, the Coulomb 
repulsion is at least U£~ l , but since there is no minimal separation, we have no direct 
handle on the possible attraction due to the field. We need a lemma, which we will 

(2) 

prove in Step 4. It concerns Eq (a), the energy of the bipolaron with no Coulomb 
repulsion, i.e., {7 = 0; 

E {2) (a) > 2E W (a) - 7a 2 /3 for all a . (3.2) 
Assuming this, the total energy in the k = region exceeds 2E^ (a) provided 

Ur 1 > 7a 2 /3 + c x t 2 , 

that is, no binding occurs if 

Ua- 1 > 7£a/3 + c x /{lot) . (3.3) 

Setting the right sides of (13. ip and (13.31) equal leads to the choice £ = c^/a and to 
absence of binding if U > Ca, as asserted. 

Step 4- The universal lower bound (13.21) : In this step, {7 = 0. We first note that 

E^(2a) > 2E {1) (a)-4a 2 /3. 

This follows from the lower bound E^(a) > —a — a 2 /3 in [23] and the upper bound 
E^(a) < —a in [HI [III 15] , stated above. So (13. 2p will follow if we can prove that 

E {2 \a) > E^(2a) -a 2 . (3.4) 

For this purpose we go back to the functional integral (I2.2p and use Schwarz's inequal- 
ity (e a+b ) < (e 2a ) 1 / 2 (e 26 ) 1//2 , where (■) now denotes expectation with respect to Wiener 
measure. We choose a to be the sum of the two terms i = j = 1 and i = j = 2 in 
(11.31) . and b to be the mixed terms i ^ j. Since (e 2 ") 1 / 2 ~ 6 -teW (2a) f or j ar g e 
inequality (13.41) will be achieved if we can show that (e 2b ) 1/2 ~ e T ° 2 . At first sight, 
the double path integral (e 2b ) looks like that for a positronium-like atom, i.e., two 
particles attracting each other through a Coulomb force with coupling constant 4a. 
The trouble is that the interaction in (11.31) is at different times, i.e., |xx(t) — x 2 (s)| _1 . 
A simple application of Jensen's inequality, however, shows that we can fix the time 
difference u = t — s and obtain the bound 



oo 



Xi(t) - x 2 (t - u)\ 
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Because of the T-periodic time translation invariance of the Wiener measure, the path 
integral is, in fact, independent of u. Hence we get the positronium-like answer as 
a bound. This completes our argument for the universal bound ( 13. 2ft and hence the 
absence of bipolaron binding for sufficiently large U/a. 



4. POLARONS IN ONE DIMENSION 

In 1976 E. P. Gross [10] wrote a seminal paper on the polaron in which he discussed 
a one-dimensional version. Even though it is not very physical, this model has been 
widely studied [HH EH [32] and we are able to prove an interesting theorem about it 
which we report here for the first time. While we have ignored the Fermi statistics up 
to now, it will play an important role in this section. 

There are N particles on the real line at xi,...,xn G R- We assume they are 
fermions, but with q spin states for each particle. The case q = N is equivalent 
to saying that Fermi statistics is irrelevant, i.e., one is dealing with boltzons. The 
Hamiltonian is as in (12.11) . except that |k| _1 is replaced by 1; the Coulomb repulsion is 
thus replaced by the delta function, and the corresponding pair potential is replaced 
by U J2i<j d( x i ~ x i)- in °ne dimension the delta function is a perfectly good potential 
of a Schrodinger operator. 

In this case we can also consider the Pekar approximation, whereby only variational 
functions of the form \1/ = ip{z\, . . . , zjv) ■ $ are allowed. Here $ is a vector in Fock 
space and Zj = (xj, af) G R x {1, . . . , q} is a space-spin coordinate for an electron. 

After minimizing the energy with respect to $, one obtains the iV-particle Pekar- 
Tomasevich functional (with spin) 

N 

I \Vi^\ 2 dZ + uY] [ 5(xi- Xj)\i){Z)\ 2 dZ - a [ p^{xf dx . (4.1) 

i=l J i<j J 

Here J dZ = a N In" ' Jm ' ' ' dx^, and the density is defined by 

p^(x)=N V" / \ip(x 1 ,a 1 ,...,x N _ 1 ,a N _ 1 ,x,(T N )\ 2 dx 1 ---dx N - 1 . 

We denote by E < ^ I \a 1 q) the infimum of ( 14.11) over all antisymmetric g-state functions 
ijj with J \ip\ 2 dZ = 1. This minimization problem also makes sense for U = oo, where 
any finite energy wave function if)(zi, . . . , Zn) must vanish if Xj = Xj for any i =fi j. We 
shall prove two facts about this minimization problem. 

Theorem 4.1. If U = and N/q is an integer, then 

E^\a,q) = (N/q)E^(a,q). 

If U = oo, then 

EW{a,q) = NEi 1 \a,l). 
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The field can cause multi-particle binding. A corollary of our first result is that, in 
the absence of repulsion, the energy per particle in the g-on state is at least as low as 
in any other state. That is, for any N (not necessarily divisible by q) 

N- 1 Ei N \a,q)>q- 1 E { *\a,q). 

To see this, consider the particle number M = Nq and apply Theorem 14.11 to this 
case. As a variational candidate for Eq (a, q) we can take q lowest energy iV-particle 
states infinitely separated from each other. Then we have Eq (a,q) < qE$ (a,q). 
On the other hand, by Theorem 4.1, E < ^ /I \a,q) = NE^\a,q), and this proves our 
assertion. Thus the g-on plays a similar role to that of nickel-62 in the curve of nuclear 
binding energies. 

When U = oo, the situation is even more dramatic; there is no binding of any kind. 

One may say that in one-dimension antisymmetry trumps the attraction caused by 
the field. (This is not true in higher dimensions.) Presumably there are finite critical 
values of U such that p-ons break apart into r-ons with p > r > 1, but we are not 
able to prove this. There should also be a finite critical value of U above which there 
is no binding of any kind. 

We now turn to the proof of Theorem 14. 1L We observe that the energy of a q-on 
can be computed explicitly, as follows. 

Lemma 4.2. If N = q, then 

E^ 9 \a,q) = -a 2 q 3 /12 . 

Proof of Lemma \4-^\ Whatever p^ might be, the minimum kinetic energy is realized 
by the product function ip(xi, . . . , x q ) = (p(xx) • • ■ (p(x q ), where f{x) = \/ p^{x)/q. 
Thus [12] 

V f \V^\ 2 dZ> [ \Vy/p-~\ 2 dx. (4.2) 

Because there are q spin states, there is an antisymmetric spin function of q variables 
with which this product function can be multiplied to yield a valid antisymmetric 
space-spin function. Equality in ( 14. 2ft is then achievable. 
To evaluate E^\a 7 q) we have to find 

E { q \a,q) = inf (q\ V '\ 2 - aq 2 \^\ 4 ) dx : \\<p\\ 2 = l| . 

The function (p(x) = (ag/4) 1 / 2 (l/ cosh(agx/2)) is easily seen to be a solution to the 
corresponding Euler-Lagrange equation and, indeed, one can prove that it is the unique 
solution of the above minimization problem (up to translations and a complex phase) 
[13j [22] . This leads to the desired expression for the energy. □ 

We need a slightly unorthodox version of a Lieb-Thirring inequality, which has been 
used before in [T7] : 
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Lemma 4.3. Assume that N/q is an integer and letip be a normalized, antisymmetric 
q -state function. Then 



N n o / /> \ 2 

\2 



Ptp(x) dx J . (4.3) 



Proof of Lemma \4-3[ Let V = —W be a negative potential in L 2 (R) and denote the 



eigenvalues of the one- dimensional Schrodinger operator —^ — W by Ai < A 2 < . . .. If 
there is only a finite number M of negative eigenvalues, we set \m+i = = . . . = 0. 

By the variational principle (see, e.g., [T8(, Thm. 12.5]) we have 



N „ N/q 

J2 / (iv^r-w^Mz)i 2 ) dz> q j2^ 

i=l J j=l 



By Holder's inequality and the sharp Lieb-Thirring inequality [22] for 3/2- moments 
of the eigenvalues, we find 



3=1 \ 3 

To summarize, we have shown that 

JV 



\ 2 /3 



J2 \ \^M 2 dZ> [ W(x)p4x)dx- N 1/3 q 2/3 [ W(x) 2 dx] 

for any < W G L 2 (W). By choosing W = cp^ and optimizing over the constant c, 
we obtain the desired bound f)4.3p . □ 



Proof of Theorem \4-l\ The case U = 0. We substitute the bound (14. 3 p into the expres- 



sion (14. ip for the energy. This lower bound only depends on the unknown quantity 
I = J Pip{x) 2 dx. By minimizing this expression with respect to / we arrive at the 
lower bound E^ N \a,q) > —a 2 Nq 2 /12. According to Lemma [4.21 this coincides with 
NE^(a,q)/q. 

To conclude the proof, we need an upper bound of the same kind. This is easily 
done by noting that we can make a state of N/q widely separated q-ons. In the limit 
that the separation goes to infinity we obtain the upper bound of N/ q times the energy 
of a single q-on. 

The case U — oo. This case is easy in view of what we just proved. The electrons, 
regardless of their spin, cannot get past each other, i.e., the iV-particle wave function 
vanishes whenever = Xj for some i ^ j. The configuration space is thus decomposed 
into a union of simplices of which S = {xi < x 2 < . . . < xjy} is representative. 
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Given a normalized g-state wave function ip we define a normalized, antisymmetric 
1-state wave function ip as follows: For x G S we set 



CT1,...,(TJV 



and we extend ?/> antisymmetrically to the other simplices. A similar construction is 
used in [19]. The crucial point is that if ip has finite kinetic energy and vanishes on the 
boundaries of the simplices, then if> has finite kinetic energy as well and vanishes on 
the boundaries of the simplices. Moreover, by the convexity inequality for gradients 
Thm. 7.8] we have 

N N 



[ \^M 2 dx<Y, [ \V^\ 2 dZ. 

8=1 1=1 ^ 



On the other hand, p^ = p^, and therefore the total energy of ip is bounded from 
below by that of ip. Note that these two energies coincide if the original ip was an 
antisymmetric function of space times a symmetric function of spin. To summarize, 
we have shown that E<£ (a, q) = E<£ (a, 1). Note that in the q — 1 case, the repulsion 
energy vanishes because the antisymmetry forces the wave function to vanish on the 
boundaries of the simplices. Thus E<£ (a, 1) = Eq(oc, 1), and the conclusion follows 
from the first part of the theorem. □ 
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